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The effect of impurities on the grain boundary melting of ice is investigated through an extension 
of Derjaguin-Landau-Verwey-Overbeek theory, in which we include retarded potential effects in a 
calculation of the full frequency dependent van der Waals and Coulombic interactions within a grain 
boundary. At high dopant concentrations the classical solutal effect dominates the melting behavior. 
However, depending on the amount of impurity and the surface charge density, as temperature 
decreases, the attractive tail of the dispersion force interaction begins to compete effectively with 
the repulsive screened Coulomb interaction. This leads to a film-thickness/temperature curve that 
changes depending on the relative strengths of these interactions and exhibits a decrease in the film 
thickness with increasing impurity level. More striking is the fact that at very large film thicknesses, 
the repulsive Coulomb interaction can be effectively screened leading to an abrupt reduction to zero 
film thickness. 
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I. INTRODUCTION 



The original microscopic theories of melting focused on 
the criterion for the homogeneous breakdown of a har- 
monic solid lattice Elaborations of this perspective 
involve the thermal activation of defects or the instability 
of the lattice I3| but a theory of melting must treat bulk 
solid-liquid coexistence and recognize that solid matter 
in laboratory and natural conditions is finite and hence 
its bulk is disjoined from the surroundings by surfaces 
0, Hence, the melting of a finite crystal begins at 
one of its free surfaces at temperatures below the bulk 
transition IJ. This interfacial premelting occurs at the 
surfaces of solid rare gases, quantum solids, metals, semi- 
conductors and molecular solids and is characterized by 
the appearance of an interfacial thin film of liquid that 
grows in thickness as the bulk melting temperature, T^, 
is approached from below. The relationship between the 
film thickness and the temperature depends on the na- 
ture of the interactions in the system. When interfacial 
premelting occurs at vapor surfaces it is referred to as 
surface melting and when it occurs at the interface be- 
tween a solid and a chemically inert substrate it is called 
interfacial melting. When films at such solid surfaces 
diverge at the bulk transition the melting is complete 
but where retarded potential effects intervene and atten- 
tuate the intermolecular wetting forces the film growth 
may be blocked and thereby be finite at the bulk transi- 
tion. This latter circumstance, in which the behavior is 
discontinuous, is referred to as incomplete melting. The 
importance of grain boundary melting is great because of 
its potentially central influence on the sintering, coarsen- 
ing, transport behavior and many other bulk properties 
in ostensibly all materials. This fact has certainly not 
gone unrecognized, but nonetheless, the great difficulty 
of directly accessing a grain boundary in thermodynamic 
equilibrium has resulted in a dearth of experimental tests. 

It is observed that, at temperatures below T,„, poly- 



crystalline matter is threaded by the liquid phase driven 
solely by the impurity and curvature depressions of the 
freezing point (e.g., 3] for ice and for metals). For ex- 
ample, the unfrozen water is found in microscopic chan- 
nels, with a scale of 10-100 /zm, called veins where three 
grains abut and at nodes separating four grains rather 
like the plateau borders in foams |3j. Hence, the equi- 
librium structure of polycrystalline ice is characterized 
by a connected network of water. Nye and Frank Q pre- 
dicted the geometry of the network under the assumption 
that the interfacial energies are independent of crystal- 
lographic orientation. In exact analogy with the clas- 
sical force balance used to determine the contact angle 
of a partially wetting fluid on a substrate or the menis- 
cus in capillary rise (e.g., 0), their result considers the 
force balance at the trijunction where three grains come 
together. This provides an expression for the dihedral 
angle, 29o, into which water intrudes, in terms of the 
simple ratio of the grain boundary, j^s , and solid-liquid, 
jsii interfacial energies: 



2cos6lo = — , 



(1) 



thereby determining the shape and cross section of a vein. 
Therefore, the concept of interfacial premelting and the 
vein-node network are related through grain boundary 
melting. Thus, during complete grain boundary melting 
the dihedral angle vanishes and the single grain bound- 
ary is replaced by two interfaces; "/ss — > '2jsi ■ This 
process is equivalent to the growth of a liquid film at the 
grain boundary and hence connects the dihedral angle, 
the vein/node/trij unction network and grain boundary 
melting . 

Because most materials exist in a polycrystalline state, 
their continuum properties, ranging from the mechanical 
and thermal properties of glaciers [^,1^ to the reduc- 
tion in the critical current density in high temperature 
superconductors 0, are to a large extent controlled by 
the character of grain boundaries. However, although 
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direct experimental probes of surface and interfacially 
melted films in thermodynamic equilibrium are common, 
experimental studies of grain boundary melting are cen- 
tered on dihedral angle measurements. The disparity in 
the relevant length scales of grain boundary films and 
the dihedral angle is such that while optical microscopy 
may be sufficient to determine the latter, the former re- 
quires a more highly resolved probe. Grain boundary 
melting poses serious experimental challenges because of 
the difficulty of direct access to the equilibrium interface 
of a bicrystal in a manner free of compromises associated 
with the proximity of surfaces associated with appara- 
tus. Although computer simulations and theoiy support 
the notion of disorder at a grain boundary ex- 
perimentally, it is often stimulated using impurities and 
then quenching a system before probing the boundary 
(e.g., |l2j |. [l^'). Experiments in aluminum using electron 
microscopy have shown that as the temperature rises, 
the boundary structure remains epitaxial until T ~ 
when the grain boundary has the signature of bulk liquid 
[l^. Dihedral angle experiments on bismuth bicrystals 
placed in contact with the melt were performed for var- 
ious crystallographic mismatches by Glicksman and col- 
leagues (See jl3 and refs therein) . They interpreted their 
data to indicate a discontinuous grain boundary melting 
transition as a function of grain mismatch. Relative to 
the abundant work on surface and intcrfacial melting, the 
present experimental edifice for grain boundary melting 
in all materials is lacking. From the standpoint of ma- 
terials science, the advantage of studying ice is that the 
melting temperature is easily accessible in the laboratory 
without substantial cryogenics or heaters. Moreover, ice 
is transparent and birefringent and hence amenable to 
optical probing. Therefore, the dual importance of grain 
boundary melting in materials in general and in ice in 
particular make a compelling case for systematic mea- 
surements in ice. Ice provides the ideal test bed for such 
measurements and its propitious geophysical importance 
is a prime motivator jig . The theory described here is 
intended to form a rostrum for experimental studies. 



II. THEORY 

A. Background 

Complete interfacial melting is determined by the com- 
petition between bulk and surface free energies and re- 
quires that the total excess surface free energy per unit 
area, -Ftotai('^)7 be a positive monotonically decreasing 
function of the film thickness with a global minimum 
at infinite film thickness. In grain boundary melting 
^totai(rf)=7ss(rf)- It is conceptually useful to consider 
the thickness dependent contributions to 7ss(<i) to arise 
from ostensibly short- and long-ranged intermolecular 
interactions; wetting forces. Hence, we write the excess 
surface free energy as ^ss{d) = 2"fse + Fshoitid) + Fiongid), 
where as above is the interfacial free energy per unit 



area of the ice-water interfaces, with implicit reference 
to the crystallographic orientation present at an inter- 
face. Hence, by definition, at large enough distances, 
the long range interactions always dominate over shorter 
range interactions. Therefore, if interactions underlying 
Fiong{d) are attractive, that is they are represented by 
a negative monotonically increasing function of d, then 
7ss(rf) can never have a global minimum at d — oo. 
However, if a short range interaction Fshort(rf) favored 
grain boundary melting, it may prevail over Fiong('i) out 
to large film thicknesses; |Fshort(rf')l ^ |^iong(rf')l 
< d < d' . The magnitude of the global minimum 
in 7ss(rf) may then be quite small but will nonetheless 
occur at a large (but finite) value of c? = dmim that 
is, 7ss(rfmin) = min [7^3 (rf)] = r. If this were the case, 
most physical observations of the system would be virtu- 
ally indistinguishable from those expected for complete 
melting. A complete calculation of jssid) allows one to 
estimate both the dihedral angle and the temperature 
at which the film thickness would saturate; each as a 
function of the thickness dmin, The dihedral angle 9o is 
given by cos 6*0 = 1-1- ^^'^ hence, in the pure case, 
we estimate the temperature at which the film thickness 
saturates as given by T — Tm « ^Tg^T"' "^here ps and 
Qm are the density and heat of fusion of ice. As we shall 
see, in the general case there is no strict separation of the 
thickness dependence of ^ss (d) because Fshort {d) can be 
quite long range, which makes the situation unique. 

Our theoretical approach employs the full frequency 
dependent dispersion (van der Waals) force contribu- 
tions to the long ranged interaction (e.g., 0, 0|) and 
Poisson-Boltzmann theory for screened Coulomb interac- 
tions. The same basic physical phenomena form the ba- 
sis of Deriagu in-Landau-Verwey-Overbeek (DLVO) the- 
ory 19, 20,^1]. The principal difference between this ap- 
proach and the modified DLVO theory used previously to 
study surface and interfacial melting ^2 is that the van 
der Waals interactions are calculated rigorously rather 
than using the phenomenological treatment of relevance 
in that case in which dispersion forces could be treated 
in the nonretarded limit. 

There are two motivations for taking a more com- 
plex approach. First, when the only interactions in 
the system are non-retarded dispersion forces, or van 
der Waals forces, to the total excess surface free en- 
ergy per unit area is most commonly denoted by 
Fvdw{d) — — TT d^, where Ah is the Hamaker 

constant (e.g., 0). For identical substrates (e.g., 
ice/ water /ice) separated by distance d the Hamaker con- 
stant is positive, producing an attraction and conse- 
quently, as has been pointed out by other authors p^l23l |. 
grain boundary melting under dispersion forces must be 
incomplete. If however the substrates are separated by 
an electrolyte solution, they may be held at bay by repul- 
sive screened Coulomb interactions. Second, for dissimi- 
lar materials, the van der Waals contribution can be both 
attractive and repulsive and one can observe oscillations 
in the force versus distance/film thickness curve leading 
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to the possibiUty of the system bein g trapp ed in a local 
rather than global minimum [13, [H, |2J, |23 ■ Such is the 
case of the surface and interfacial melting of ice 0, 0| . 
We do not expect this particular behavior in the case 
of the grain boundary under the influence of dispersion 
forces alone. However, when combined with repulsive 
screened Coulomb interactions, the effective interfacial 
free energy does indeed exhibit behavior that is analo- 
gous to the dispersion force contribution to the interfa- 
cial melting of ice. For example, as discussed in general 
terms above, depending on the electrolyte concentration, 
we can observe a global minimum at infinite, finite, or 
zero film thickness, but with a maximum at intermediate 
film thickness and a positive monotonically decreasing 
free energy at large film thickness. These phenomena, in 
part associated with retardation, thereby heightening our 
appreciation of the importance of a quantitative treat- 
ment of the long ranged forces in the system requiring a 
more complete calculation of F^j^{d). 

The calculation of frequency dependent dispersion 
force contributions to grain boundary melting are de- 
scribed in section FlI CI but it is prudent to first remind 
ourselves of the general understanding of ice surfaces un- 
der the influence of such interactions. Generally speak- 
ing, when dispersion forces dominate, the wetting of any 
ice surface by water at temperatures below Tm will be fa- 
cilitated when the polarizability of the water lies between 
that of the ice and the other material be it a gaseous 
phase, a chemically inert solid or ice. A particular nov- 
elty of the system, first pointed out by Elbaum and Schick 
[l7| in their study of the surface melting of ice, is that the 
appropriately-transformed polarizability of ice is greater 
than that of water at frequencies higher than approxi- 
mately 2 X 10^^ rad s^ ^, whereas it is smaller at lower 
frequencies. Thus, so long as the surface melted layer of 
water is thin, the polarizabilities at all frequencies con- 
tribute additively to F^dw(d)i whereas upon thickening 
retardation-the attenuation of the interaction due to the 
finite speed of light-reduces the high frequency contribu- 
tions and favors those in which the polarizability of water 
dominates over that of the ice. Hence, with regard to the 
dispersion force contribution to the premelting of an ice 
surface, we expect the process to be complete when the 
third phase in the system possesses a low frequency po- 
larizability greater than that of water, and unattainable 
if the low frequency polarizability is less than that of wa- 
ter. Many, if not most, solids have dielectric properties 
that lead to complete interfacial premelting whereas the 
vapor phase of ice does not. 

The great sensitivity of surface meltin g ex periments to 
the presence of impurities led one of us |22l| to introduce 
electrolytes into a model that did not include retarda- 
tion. Here, as in that case, we consider the ions that 
underlie the repulsive screened Coulomb interactions in 
the canonical ensemble and the entire ice/solution system 
is treated in the grand canonical ensemble Hence, 
because the ions are insoluble in the solid phase and the 
Debye length characterizing the range of the interactions 



is inversely proportional to the square root of the ion 
density, then the range and strength of this contribution 
to grain boundary melting is sensitive to the impurity 
level. Therefore, it is required that we properly capture 
the competition between repulsive screened Coulomb in- 
teractions and attractive dispersion forces at very long 
range, which is done through the use of Debye-Hiickel 
theory in the appropriate ensemble and the complete the- 
ory of Dyzaloshinskii, Lifshitz and Pitaevskii, We 
now describe the various contributions to the excess sur- 
face energy in turn, the total free energy of the system, 
and the associated predictions. 

B. Ionic Force Effects 

The screening of surface charge is described by the 
Poisson-Boltzmann (T*B) equation for the electrostatic 
potential ip created by the distribution of ions of num- 
ber density n{z) throughout the film. For a monovalent 
electrolyte, and a surface potential ips less than 25 mV, 
the Poisson-Boltzmann equation can be linearized in the 
Debye-Hiickel limit which yields 

i^iz) = ?^se-"^ where = ( ^ (2) 

is the Debye length, which captures the characteristic fall 
off of the ion field in the direction z normal to the surface. 
The bulk ion density is n;,, e is the dielectric constant of 
the film, Eo is the free space permittivity, e is the elemen- 
tary charge, and kbT is the usual thermal energy. The 
repulsive force between two charged surfaces originates 
in the restriction of the configurational entropy of the 
ions as the surfaces are brought closer. The resulting ex- 
cess interfacial free energy per unit area across a film of 
thickness d is written 

FoHid) « ^e-'^ , (3) 

K€€o 

where Qs is the surface charge density, discussed in more 
detail in section IH Dl In the classical context of DLVO 
theory, the Debye length is a constant, and equation Q 
can be considered in the same sense as i^short(<^) above. 
At this stage we simply note that it describes the electro- 
static contribution to total excess interfacial free energy 
of relevance to the grain boundary melting of ice in the 
presence of electrolytes. 

We point out several important modifications to PB 
theory beyond the obvious simplifications associated with 
the linearization, or weak coupling limit (e.g., 28]), that 
led to equation |5J. First, as a mean field theory, ion- 
ion correlations and discrete ion-solvent interactions are 
not captured and the solvent is treated as a continuous 
medium 12^ . Second, there are dispersion forces between 
ions and the surfaces that bound them, which can be of 
much longer range than Debye-Hiickel theory would pre- 
dict although not of significant magnitude at a range be- 
yond about 3 nm where retardation becomes important 
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[2l|. Such ion dispersion forces can influence the static 
van der Waals interactions in the system significantly for 
high dopant levels or multivalent ions 30]. Third, at 
molecular length scales steric, or finite ion size, effects are 
important (31.] . Although all of these phenomena may be 
of importance in a limited regime of the predictions pre- 
sented here, we do not expect them to be responsible for 
a qualitative difference the principal features described. 
Moreover, other equally distinct effects can also play a 
role. For example, in the case of the ion dispersion phe- 
nomena, the fact that the chemical potential of the ions 
in ice is effectively zero results in the film concentrat- 
ing or diluting as a function of the departure from bulk 
coexistence leads to a film thickness dependence of the 
Debye-length. Hence, we find a very long range repulsive 
interaction with an effectively algebraic structure. As for 
steric infiuences, at a grain boundary other proximity ef- 
fects associated with the crystallographic mismatch will 
be come equally important at short range. A theory that 
systematically includes all such effects on equal footing is 
still unwarranted by the state of the experimental land- 
scape, but we are compelled to note for the reader that 
we are acutely aware of all such complications and lim- 
itations enumerated in the detailed studies cited in this 
paragraph and in books on the topic jl^], 1^3- 



C. Dispersion Force Effects 

The frequency dependent dispersion force, or van der 
Waals, contribution to the free energy of a surface con- 



sisting of a pure liquid water layer of thickness d between 
two slabs of bulk ice can be obtained from the theory of 
Dzyaloshinskii, Lifshitz, and Pitaevskii [23] (DLP). The 
result of this theory is an integral expression for the in- 
terfacial free energy per unit area, F^dwid), in terms of 
the frequency-dependent dielectric polarizabilities of the 
ice(s) and the water (£). In the sense described in section 
III Al were grain boundary melting under dispersion forces 
alone to be considered, then 7ss (d) = 27^^ -I- Fiong(rf) = 
27s£ + ^•udw (d) ■ Were complete grain boundary melting 
to be expected, and it is not, it would be indicated by a 
global minimum of Fy^w (d) at d — > 00 so that — ^Isi 
at bulk coexistence. Incomplete grain boundary melt- 
ing under dispersion forces is indicated by a minimum at 
finite rf, with F^dwid) negative there. 

As described in III Al retarded potential effects are 
known to influence the qualitative nature of melting in 
ice. For example, the surface melting of ice has been 
observed to be incomplete due to retarded potential ef- 
fects, whereas slight impurity doping leads to complete 
surface melting |17l l2a ] . Because the neglect of retarda- 
tion can lead to qualitatively and quantitatively different 
predictions we append a brief extended discussion of the 
matter. 

In order to make this development reasonably self con- 
tained we present the DLP expression for F^dwiji) de- 
spite the fact that is has appeared in a variety of forms 
Ullll I2II23, and with sympathy to the girth of the 
literature we attempt the utmost brevity. The frequency 
dependent dielectric response of the ice/water/ice system 
is described by 
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(5) 



and the material (s, £) dielectric functions, corresponding 
to ice (e^) and water (e^) are evaluated at the sequence 
of imaginary frequencies — i(2'KkT /h)n. The prime 
on the sum, indicates that the n = term is weighted by 
1/2. The lower limit of integration is r„ = 2(i(e£)^/^^„/c, 
and fc, h and c have their usual meaning. The dielectric 
function required in the integral, e(jO, is the analytic 
continuation of the material dielectric function e{uj) to 
imaginary frequencies. Lacking complete spectra for ice 
and water, as can be obtained in high melting temper- 
ature materials (e.g., [I3j ]) we must generate the func- 
tion by fitting the dielectric response of the material to 



a damped-oscillator model of the form, 

ii 



(6) 



where e^, fj and gj are fitting parameters |33j]. Each 
term in the sum corresponds to an absorption band of 
frequency, width and oscillator strength, e^, gj and 
respectively. Substitution of for uj gives e(jC)i ^ "weW 
behaved, monotonically decreasing, real function of ^. 

We note that equation assumes implicitly that is 
isotropic. Although the effect of the crystal orientation 
has been studied theoretically |23], no such polarization 
data was available for the present calculation. It has not 
been measured over the full frequency range, but judg- 
ing from known values at optical frequencies |35l ]. it is 
estimated to be small, and hence we neglect the small 
orientational dependence of the polarizabilities. There- 
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fore, we treat the ice crystals as continuum dielectrics, 
with the surface crystal structure embodied by the sur- 
face charge density as described in section HTbI 

D. Total Free Energy 

The stable existence of an impure grain boundary film 
at temperatures below is investigated by minimiza- 
tion of the total free energy of the system. The distinc- 
tion between the canonical DLVO setting, say in colloid 
science, and surface, interfacial or grain boundary melt- 
ing is that in the former case the Debye length is constant 
for all values of the film thickness whereas this is not the 
circumstance in the latter case We imagine doping 
a pure film with a monovalent electrolyte such as NaCl. 
Depending on the nature (low angle or twist) and mag- 
nitude of the mismatch between the crystals flanking the 
grain boundary, the solid-liquid interface will have a par- 
ticular surface charge density Qs- We realize that each 
solid-liquid interface can, in principle, have a different 
value of qs, and we embody the surface crystal structure 
through the surface charge density, but for the sake of 
simplicity and pedagogy we assume that both are equal. 
The surface charge is screened by the ions in the film, 
creating a repulsive interaction which competes with at- 
tractive long ranged van der Waals interactions and this 
embodies the excess surface energy. 

We deposit Ni moles per unit area of a monovalent 
electrolyte in the grain boundary. The chemic al p oten- 
tial of the electrolyte in ice is essentially zero |33, and 
because we are well below the solubility limit the ions 
remain in the film with a concentration that is inversely 
proportional to film thickness. When the temperature 
decreases the film thins, the ion concentration increases 
and hence the Debye length decreases (see equation (|2l). 
At very low film thicknesses and very high ion concen- 
trations the limitations of PB theory, discussed in section 
111 Bl come to the fore. Within these approximations the 
behavior of the total free energy of the system determines 
the nature of grain boundary melting. For this we com- 
bine bulk and surface terms as 

Gt{T, P, d, N,) = pi^iiiT, P)d + ^i,{T, P)N, 

+RTN, en^+ M.(r, P)Ns + 7s. (d), (7) 
Pea 

where jss{d) = 2jsi+Fi:iH{d)+Fydw{d) is the total effec- 
tive interfacial free energy combining equations ^ and 
The molar density, chemical potential per mole and 
the number of moles per unit area of the solvent are pi, 
pi and Ne—pid respectively and the chemical potential 
of the impurity, the solid and the number of moles of the 
latter are /z^, ps and Ng. We write the mixing entropy 
term, RTNi In where R is the gas constant, in the 
dilute limit of ideal solution theory although we could 
simply replace the molar ratio by the activity coefhcient 
for the impurity in the solvent. These complications do 
not change the principal regimes exhibited by the system. 



The total free energy is minimized with respect to film 
thickness at fixed temperature and dopant concentration 
and the integral in Eq. and the dielectric fits for 
ice and water compiled from reference [36|, were calcu- 
lated in the same manner as was done previously p^.p7||. 
The Gauss-Legendre quadrature integration scheme was 
tested against these calculations for the case of surface 
melting in the absence of impurities. 



III. RESULTS 

A representation of what we see overall as we vary 
the parameters is shown in the figure. In a given cal- 
culation, we ascribe a surface charge density to both 
solid/film interfaces, with, for example qg — 0.01 C 
corresponding to approximately four elementary charges 
per 1000 molecular sites. This particular case is rather 
high from the perspective of counting missing bonds at 
the interface or ionization defects, but conservative from 
the perspective of charge adsorption. It is useful to be- 
gin in the upper left hand corner of figure 1(a) at high 
temperature. Consider the most dramatic curve first; 
the solid line which has the lowest dopant level of 6 x 
10^^/iMm^'^. Due to the fact that in this case the De- 
bye length is extremely long and the repulsive Coulomb 
interaction is longer ranged than the attractive van der 
Waals interaction, a very thick film of liquid exists at 
the ice grain boundary. As the temperature is lowered, 
the efficiency with which ice rejects electrolytes causes 
the film to be enriched with ions, and hence the Debye 
length, and the range of the Coulomb interaction, de- 
creases until such point that the attractive tail of the 
van der Waals force begins to come into its own. As the 
temperature decreases a point is reached where abruptly 
the van der Waals interaction dominates and drives a 
discontinous transition to zero film thickness. The tran- 
sition moves to higher temperatures as the surface charge 
density and dopant levels decrease thereby showing the 
expected behavior for melting driven solely by dispersion 
forces, that is, no grain boundary melting. An increase in 
the initial dopant level reduces the Debye length and the 
strength and range of the screened Coulomb interaction. 
Hence, the plateau where the van der Waals and screened 
Coulomb interactions compete shifts to lower film thick- 
nesses. If the temperature and/or the dopant level is suf- 
ficiently high, the bulk free energy dominates the total 
free energy and, as found previously a surface ver- 
sion of Raoult's law is seen; d oc (T^ — T)~^. Finally, the 
largest dopant levels lead to the bulk free energy domi- 
nating through the entire range, and we predict a measur- 
able film thickness in an easily accessible and controllable 
experimental range of temperature from to - 10°C. Al- 
though this latter behavior was observed in studies of 
interfacial melting using atomic force microscopy |37j . 
there exist no observations on ice grain boundaries. We 
emphasize that at fixed temperature the film thickness 
is not a monotonically increasing function of ion concen- 
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tration until a threshold level is surpassed. Thus, there 
is a transition from melting dominated by dispersion and 
ionic force interactions to that dominated by the more 
commonly considered colligative effects. An exploration 
of variations in dopant levels and surface charge density 
reveals the same qualitative behavior. 
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FIG. 1: Film thickness (A) versus undercooling (Kelvin) for 
various monovalent electrolyte (NaCl) concentrations at an 
ice grain boundary with a solid/film surface charge density Qs 
of (a) 0.01 C (b) 0.001 C and (c) 0.0001 C m'^ In 

(a)-(c) Ni = 6x 10"^ (solid line), 6 x 10"* (long dashed line), 
6 X 10~^ (medium dashed line), 1.8 x 10"^ (small dashed 
line) and 6 (dotted line) ^Mm~^. At the lower dopant levels 
where the Debye length is long we observe an abrupt or dis- 
continuous onset of grain boundary melting. The transition 
moves to higher temperatures as the surface charge density 
and dopant levels decrease thereby showing the expected be- 
havior for melting driven solely by dispersion forces. Note, for 
example in (a), that at Tm —T — 0.01 K, the film thickness is 
not a monotonically increasing function of ion concentration. 
Indeed, the film thickness decreases with increasing concen- 
tration until a threshold level is surpassed. Hence, there is 
a transition from melting dominated by dispersion and ionic 
force interactions to that dominated by the more commonly 
considered colligative effects. In (b) and (c) we see that for 
low surface-charge densities the Coulomb interaction is weak, 
resulting in a much less-pronounced plateau, or complete ab- 
sence thereof, except for the lowest dopant concentrations. 



IV. DISCUSSION 

One effect not considered in the current theoretical 
framework that warrants discussion is the possible depen- 
dence of jsi on impurities adsorbed at the interface. The 
low solubility of electrolytes in ice leads to an enhanced 
interfacial adsorption and will modify the solid/film in- 
terfacial free energy, jsi, in a manner that depends on, 
among other things, the magnitude of the surface charge 
density and the bulk electrolyte concentration. A de- 
crease in "fsi will enhance grain boundary melting rela- 
tive to the pure case, were there not a similar decrease 
in the dry grain boundary energy. Although it is well 
known that an image charge effect leads to an increase 
in the liquid/ vapor interfacial free energy with electrolyte 
concentration, there is no quantitative measure of the 
phenomenon at the solid/liquid interface wherein the in- 
trinsic ionization defects in the solid phase are available 
to screen the adsorbed ions. Indeed, the great sensitivity 
of the surface melting of ice to impurities, in part respon- 
sible for the transition from incomplete to complete melt- 
ing p^ . is quite strikingly demonstrated in the calcula- 
tions shown here. However, deconvolution of the precise 
dependencies of the constituents of the interfacial ener- 
gies on Ni, the coupling to the ionic forces in the solution 
within the PB framework, and the additional influences 
of ion-ion correlations, ion-solvent and ion-dispersion in- 
teractions, and steric effects described above, is beyond 
the scope of experimental evidence. 

As described above, there are two principal ways in 
which grain boundary melting can be observed; scrutiny 
of the dihedral angle subtended by a grain boundary 
groove, or direct probing of the liquid film thickness at 
the boundary as a function of temperature. In principal, 
the dramatic nature of our predictions and their ther- 
modynamic accessibility renders them amenable to an 
experimental search. In the experimental scenarios that 
we envisage, the grain boundary of a bicrystal is in con- 
tact with the bulk solution at a classical grain boundary 
groove and hence, it is theoretically possible to examine 
both the dihedral angle and the grain boundary thick- 
ness. However, the ensemble operative in our predictions 
leads to a particularly impurity sensitive film thickness 
dependence of the Debye length; Equation ^ leads to 
K = c y'Ni/d, where c = 7.237 x 10^ m^/^ ^o\-i/2 jg ^ 
constant. If however, contact of the grain boundary film 
with a bulk reservoir of ions allows compositional equi- 
libration between the film and the bulk, then the Debye 
length will not be a function of the film thickness and 
some qualitative aspects of the predictions will differ. 
We expect compositional equilibration for the thickest 
films, but the presence of a substrate field in the case of 
the thinner films may not facilitate such homogenization. 
This question is the subject of ongoing work. 

There are different methods that probe different as- 
pects of premelted liquid water. For example, infrared 
spectroscopy averages the premelted liquid in a polycrys- 
talline sample , introduction of a force microscope tip 
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of another material provides and indirect probe |33 , |39j , 
bright X-ray reflectivity discerns short range structural 
properties and optical methods (e.g., [2^, 
tinguish between liquid water and ice. Regardless of the 
approach, our calculations indicate the preference of sys- 
tematic doping over attempts to prepare a completely 
clean system. 

In summary, we have studied the effect of impurities 
on the grain boundary melting of ice using an extension 
of Derjaguin-Landau-Verwey-Overbeek theory, in which 
we include retarded potential effects in a calculation of 



the full frequency dependent van der Waals and Coulom- 
bic interactions within a grain boundary. We find that, 
depending on the amount of impurity and the surface 
charge density, as temperature decreases, attractive dis- 
persion force interactions effectively compete with re- 
pulsive screened Coulomb interactions producing grain 
boundary melting through a discontinuous transition. At 
sufficiently high dopant concentrations the classical solu- 
tal effect dominates the melting behavior. The effect is 
within the scope of experimental accessibility. 
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APPENDIX A: RETARDATION EFFECTS 

The detailed quantitative effects of retardation men- 
tioned in section III Al are placed in this Appendix be- 
cause grain boundary melting under dispersion forces 
alone is incomplete. Retardation has been clearly de- 
scribed previously (e.g., 0,0,13,0) by considering a 
general ice/water/X interface in which X denotes a "sub- 
strate" which may be pure water vapor. In the limit that 

« « ex ~ 1, one can approximate Eq. ^ as; 

n— 

(Al) 

The term e due to retardation, acts as a high fre- 
quency cutoff to the sum which is inversely proportional 



to d. 

When the substrate is pure water vapor, ex may be 
taken equal to 1. In this case, it is clear that were 
Es — Q < at all frequencies, Fydw{d) would be a mono- 
tonically increasing function of d and the film would not 
grow. Conversely, if ~ > at all frequencies, then 
Fvdwid) would be a monotonically decreasing function 
and the film thickness would diverge as the melting tem- 
perature is approached. Because eg — ei changes sign at 
frequency ^c, the resulting melting behavior is interme- 
diate between these cases. In particular, for sufficiently 
large d, where the sum is dominated by the low frequency 
terms, surface melting is inhibited. When the substrate 
is a solid material with arbitrary dielectric properties, 
the results are more complicated. The function ex now 
depends on frequency ^, and ei{i^) — ex{iC} ni^-Y change 
sign. For all of the materials studied previously we 
observed that ex > es^i in the frequency range ^ < ^c- 
Hence we found that for d large enough («30A) to allow 
retardation to come into play, Fy^wid) will be a positive 
monotonically decreasing function of d, implying that the 
necessary condition for complete interfacial melting to 
occur is fulfilled. To determine whether there is another, 
deeper, minimum in the free energy at small film thick- 
nesses requires implemention of the full dispersion force 
calculation embodied in equation (0J. 



